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Q ■ Abstract 

We measure the topological charge and its fluctuation for the gauge configurations generated by 
the RBC and UKQCD Collaborations using 2+1 flavors of domain- wall fermions on the 16 3 x 32 
lattice (L ~ 2 fm) with length 16 in the fifth dimension at the inverse lattice spacing a -1 ~ 
1.62 GeV. From the spectral flow of the Hermitian operator H w {2 + 75iT u ,) — we obtain the 



topological charge Qt of each gauge configuration in the three ensembles with light sea quark 
masses m q a = 0.01, 0.02, and 0.03, and with the strange quark mass fixed at m s a = 0.04. From 



> . 

our result of Qt, we compute the topological susceptibilty xt = (Q?)/^> where Q is the volume of 
the lattice. In the small m q regime, our result of xt agrees with the chiral effective theory. Using 
the formula xt = ^/(m" 1 + m^ 1 + m~ l ) by Leutwyler-Smilga, we obtain the chiral condensate 
OO : S MS (2 GeV) = [259(6) (9) MeV] 3 . 
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I. INTRODUCTION 



In Quantum Chromodynamics (QCD), the topological susceptibility (xt) is the most 
crucial quantity to measure the topological charge fluctuation of the QCD vacuum, which 
plays an important role in breaking the Ua(1) symmetry. Theoretically, Xt is defined as 



Xt = I d'x (p(x)p(O)) , (1) 

where 



P(x) = ^^e^ Xa tY[F^(x)F Xa {x)], (2) 



is the topo.og.ca. charge density expressed in term of the matrix-va.ued ne.d tensor F„. 
With mild assumptions, Witten [1] and Veneziano [2] obtained a relationship between the 
topological susceptibility in the quenched approximation and the mass of rj' meson (flavor 
singlet) in unquenched QCD with Nf degenerate flavors, namely, 

X< (quenched) = 

where f n = 131 MeV, the decay constant of pion. This implies that the mass of rf is 
essentially due to the axial anomaly relating to non-trivial topological charge fluctuations, 
which can turn out to be nonzero even in the chiral limit, unlike those of the (non-singlet) 
approximate Goldstone bosons. 

Using the Chiral Perturbation Theory (ChPT), Leutwyler and Smilga [3] obtained the 
following relation in the chiral limit 

*= 1 j\ +0{™&, (N f = 2 + 1), (3) 

where m u , m^, and m s are the quark masses, and E is the chiral condensate. This implies 
that in the chiral limit (m u — > 0) the topological susceptibility is suppressed due to internal 
quark loops. Most importantly, (J3]) provides a viable way to extract E from xt in the chiral 
limit. 

From (CQ), one obtains 

(Q?) 



Qt = J d 4 xp(x), (4) 
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where Q is the volume of the system, and Qt is the topological charge (which is an integer 
for QCD). Thus, one can determine Xt by counting the number of gauge configurations for 
each topological sector. Furthermore, we can also obtain the second normalized cumulant 

1 



C4 = "fi 



(Q^>-3(Q 2 ) 2 , (5) 



which is related to the leading anomalous contribution to the rf — rf scattering amplitude 
in QCD, as well as the dependence of the vacuum energy on the vacuum angle 9. (For a 
recent review, see for example, Ref. {4] and references therein.) 

However, for lattice QCD, it is difficult to extract p(x) and Q t unambiguously from the 
gauge link variables, due to their rather strong fluctuations. 

To circumvent this difficulty, one may consider the Atiyah-Singer index theorem 

Qt = n + — ri- = index(D), (6) 

where n± is the number of zero modes of the massless Dirac operator V = 7 M (<9 M + igA^) 
with ± chirality. 

For lattice QCD with exact chiral symmetry, it is well-known that the overlap Dirac 

n q 

operator [6|, |7|] in a topologically non-trivial gauge background possesses exact zero modes 
(with definite chirality) satisfying the Atiyah-Singer index theorem. Thus we can obtain the 
topological charge from the index of the overlap Dirac operator. Writing the overlap Dirac 
operator as 

D = m (1 + 75-^1- (7) 



where H w is the standard Hermitian Wilson operator with negative mass —mo (0 < mo < 2), 
then its index is 



index (D) = Tr 



D 

75 1- 



2mr 




n+ - n_ = Q t , (8) 



where Tr denotes trace over Dirac, color and lattice spaces. 
Obviously, from (jSJ), we have 

index(D) = n + - n_ = -^Tr l^jL j = \{h-- h + ), (9) 
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where h + (hJ) is the number of positive (negative) eigenvalues of the hermitian Wilson-Dirac 
operator H w . However, one does not need to obtain all eigenvalues of H w in order to know 
how many of them are positive or negative. The idea is simple. Since H w has equal number 
of positive and negative eigenvalues for m < 0, then one can just focus on those low-lying 
(near zero) eigenmodes of H w , and see whether any of them crosses zero from positive to 
negative, or vice versa, when ttlq is scanned from zero up to the value used in the definition 
of D. From the net number of crossings, one can obtain the index of D. This is the spectral 
flow method used in Ref. to obtain the index of the overlap Dirac operator. We also 
applied the spectral flow method to obtain the index of the overlap Dirac operator and to 
determine the topological susceptibility in quenched QCD In this paper, we extend our 
previous studies to unquenched QCD. 




FIG. 1: The spectral flow of 12 lowest-lying eigenvalues of H for the gauge configuration number 
1825 in the ensemble (m q = 0.01, m s = 0.04). There are 10 net crossings from negative to positive, 
so the index is —10. 

For the conventional domain- wall fermion [lo| , its effective 4-dimensional Dirac oper- 
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at or is 



B = ™ £ p-™ 2 )j i+75 _|_j i (10) 

where 

H = H w (2 + l 5 H w )-\ (11) 

Thus, we can obtain the topological charge of a gauge field configuration from the spectral- 
flow of if as a function of m . Obviously, it is much more computationally intensive to 
project the low-lying eigenvalues of H than those of H w , due to the extra inverse operator 
(2 + l5 H w )-\ 

In this paper, we use the spectral flow of H to determine the topological charge of the 
gauge configurations ( |http: / /lattices.qcdoc.bnl.gov/[ ) generated by the RBC and UKQCD 
Collaborations using 2+1 flavors of domain- wall fermions on the 16 3 x 32 lattice (L ~ 2 fm) 



with length 16 in the fifth dimension at the inverse lattice spacing a -1 ~ 1.62(4) GeV 11]. 
There are three ensembles of gauge configurations with light sea quark masses m q a = 0.01, 
0.02 and 0.03, and with the strange quark mass fixed at m s a = 0.04. For the ensmeble 
with m q a = 0.01, we pick one configuration every 5 configurations, from configurations 
numbering from 0020 to 4015. Thus we have 800 configurations with m q a = 0.01. Similarly, 
for m q a = 0.02, we pick 809 configurations from configurations numbering from 0005 to 
4045, and for m q a = 0.03, we pick 717 configurations from configurations numbering from 
4020 to 7600. 

In Fig. [1], we plot the spectral flow of 12 lowest-lying (near zero) eigenvalues of if (mo) in 
the interval 0.8 < mo < 1.8, for the gauge configuration number 1825 in the ensemble with 
m q = 0.01 and m s = 0.04. In this case, the net crossings from negative to positive is 10, so 
the index is —10. In general, it may happen that there are some intriguing eigenvalues lying 
very close to zero (e.g., the one around m = 1.45 in Fig. [1]). Thus, with a coarse scan in m , 
it may not be able to determine whether they actually cross zero or not. These ambiguities 
can only be resolved by tracing them closely at a finer resolution in mo- Obviously, it is a 
very tedious job to determine the topological charges of 2,326 gauge configurations via the 
spectral flow of H = H w {2 + ^H w ) -1 . 
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FIG. 2: Histogram of topological charge distribution for m q a = 0.01,0.02, and 0.03 respectively. 



II. RESULTS 

In Fig. |2l we plot the histogram of topological charge distribution for m q a = 0.01, 0.02, 
and 0.03 respectively. Evidently, the probability distribution of Q t becomes more sharply 
peaked around Q t = as the light sea quark mass m q gets smaller. Our results of topological 
susceptibility Xt, and the second normalized cumulant C4, together with their ratios c^/xt, 
and c^/(2xt^) are listed in Table [H The error is estimated using the jackknife method with 
bin size of 20 configurations for m q = 0.01 and 0.02, and 13 configurations for m q = 0.03, 
with which the statistical error saturates. 

Evidently, the statistical error of the topological susceptibility is about 10%, while that 
of C4 is very large due to low statistics. Therefore, one cannot draw any conclusions from our 
result of C4, as well as from the ratio C4/ Xt- Interestingly, our result of C4/ (2%^) is consistent 
with that in Ref. [ijj], which is obtained from the plateaus (at large time separation) of the 
2-point and 4-point time-correlation functions of the flavor-singlet pseudoscalar meson rf in 
a fixed global topology with Q t = 0. 

In Fig [3j we plot our data of Xt versus m^, where the pion mass m^a and the inverse 



lattice spacing a 1 are determined by the RBC and UKQCD Collaborations 11]. 



For three flavors with m u = rrid, we may use the partial conservation of the axial current 
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m a a 


Xt 


c 4 




c 4 /(2Xt^) 


0.01 


5.51(62) x 10" 5 


-1.88(8.33) x 10~ 5 


-0.34(1.55) 


-0.0239(1110) 


0.02 


7.74(75) x 10~ 5 


-4.99(16.18) x 10" 5 


-0.64(2.15) 


-0.0317(1088) 


0.03 


1.23(11) x 10~ 4 


4.69(3.46) x 10~ 4 


3.80(2.57) 


0.1172(834) 



TABLE I: The topological susceptibility xt, the second normalized cumulant C4, and their ratios 
Ci/xt, and C4/(2xfO), versus the sea quark masses, for Nf = 2 + 1 lattice QCD with domain-wall 
fermions. 



(PCAC) relation m 2 PS = Cm q to transcribe the Leutwyler-Smilga relation ([3j) to 

Xt = - 



1 + 



2m 2 ps (rn s ) 



(12) 



where mps(m s ) is the mass of the pseudoscalar meson with strange quarks, and £' = £C . 
The data points of Xt are well fitted by the ChPT formula $F2ft with £' = 0.0039(2) [GeV 2 ]. 
The fitted curve is plotted as the solid line in Fiejl Using C = 7.4034(1986) [GeV] deter- 
mined by the RBC and UKQCD Collaborations 11], we obtain £ = 0.0288(16) [GeV 3 ]. 

In order to convert £ to that in the MS scheme, we use the renormalization factor 
Zp(2 GeV) = 0.604(18) (55) which is determined by the RBC and UKQCD Collabora- 
tions [l^Lemploying the non-perturbative renormalization technique through the RI/MOM 
scheme 13j. Then the value of £ is transcribed to 

£^(2 GeV) = [259(6) (9) MeV] 3 , 

which is in good agreement with the results extracted from xt m 2+1 flavors QCD [14] and 



2 fl 



avors QCD 



15| . as well as that obtained from the low-lying eigenvalues in the e-regime 
16j. The errors represent a combined statistical error (a -1 and -Z^ s ) and the systematic 
error respectively. Since the calculation is done at a single lattice spacing, the discretization 
error cannot be quantified reliably, but we do not expect much larger error because the 
domain- wall fermion action is free from 0(a) discretization effects. 



III. CONCLUDING REMARK 



In this paper, we have obtained the topological susceptibility Xt and the second normal- 
ized cumulant C4, in 2+1 flavors lattice QCD with domain-wall fermions. The expected sea 
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FIG. 3: The topological susceptibility xt versus for 2+1 flavors lattice QCD with domain-wall 
fermions. 

quark mass (pion mass square) dependence of Xt from ChPT is clearly observed. However, 
our statistics of ~ 800 configurations are insufficient to determine C4 with a reasonably small 
error. 

Finally, it is interesting to compare the topological susceptibility obtained from the global 
topological charge (this work) with that extracting from the plateaus (at large time separa- 
tion) of the 2-point and 4-point time-correlation functions of the flavor- singlet pseudoscalar 
meson rf with fixed topology [14]. We find that both methods give Xt i n agreement with 
the chiral effective theory. Furthermore, the chiral condensates extracted from both sets of 
Xt are also in good agreement. Next, we turn to the second normalized cumulant C4. If 
we try to measure it with the global topological charge, it would require more than 10,000 
configurations in order to pin down the statistical error less than 10%. On the other hand, 
if we measure C4 with the correlation function of topological charges in sub-volumes (with 
fixed global topology), we would expect that it only requires about 1,000 configurations in 
order to achieve a level of 10% statistical error. To conclude, it is interesting to see that 
there are more than one viable options to obtain the topological susceptibility and the higher 
normalized cumulants, in lattice QCD with exact chiral symmetry. 
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